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Electrons most often organize into Fermi-liquid 
states in which electron-electron interactions 
play an inessential role. A well known exception 
is the case of one-dimensional (ID) electron sys- 
tems (IDES). In ID the electron Fermi-surface 
consists of points, and divergences associated 
with low-energy particle-hole excitations abound 
when electron-electron interactions are described 
perturbatively. In higher space dimensions, 
the corresponding divergences occur only when 
Fermi lines or surfaces satisfy idealized nesting 
conditions. In this article we discuss electron- 
electron interactions in 2D graphene bilayer 
systems which behave in many ways as if they 
were one-dimensional, because they have Fermi 
points instead of Fermi lines and because their 
particle-hole energies have a quadratic dispersion 
which compensates for the difference between ID 
and 2D phase space. We conclude, on the basis 
of a perturbative RG calculation similar to that 
commonly employed in ID systems, that interac- 
tions in neutral graphene bilayers can drive the 
system into a strong-coupling broken symmetry 
state with layer-pseudospin ferromagnetism and 
an energy gap. 



Recent progress in the isolation of nearly perfect sin- 
gle and multilayer graphene sheet*^^ has opened up 
a new topic in two-dimensional electron systems (2DES) 
physics. There is to date little unambiguous experimen- 
tal evidence that electron-electron interactions play an 
essential role in the graphene family of 2DES's. How- 
ever, as pointed out by Min et aZ£ graphene bilayers near 
neutrality should be particularly susceptible to interac- 
tion effects because of their peculiar massive-chiral^ band 
Hamiltonian, which has an energy-splitting between va- 
lence and conduction bands that vanishes at k = and 
grows quadratically with k = \k\: 
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In Eq. (JTJ the r z s are Pauli matrices and the Greek la- 
bels refer to the two bilayer graphene sublattice sites, 
one in each layer, which do not have a neighbor in the 
opposite graphene layer. (See Fig. [TJ) The other two 
sublattice site energies are repelled from the Fermi level 
by interlayer hopping and irrelevant at low energies. It 
is frequently useful to view quantum two-level layer de- 
gree of freedom as a pseudospin. The 3 — 1 pscudospin 



chirality of bilayer graphene contrasts with the J = 1 
chiralitjii of single-layer graphene and is a consequence 
of the two-step process in which electrons hop between 
low-energy sites via the high-energy sites. The massive- 
chiral band-structure model applies at energies smaller 
than the interlayer hopping scaled 71 ~ 0.3 eV but larger 
than the trigonal- warping scaled 7371/70 ~ 0.03 eV be- 
low which direct hopping between low-energy sites plays 
an essential role. The body of this paper concerns the 
role of interactions in the massive-chiral model; we re- 
turn at the end to explain the important role played by 
trigonal warping. 




FIG. 1: Bilayer graphene lattice structure and elec- 
tronic structure. The massive chiral fermion model de- 
scribes the low-energy sites in a AB-stacked graphene bilayer, 
those atom sites (top layer B sites and bottom layer A 1 sites) 
which do not have a neighbor in the opposite layer. The con- 
duction and valence bands touch at the Brillouin-zone corner 
wavevectors, taken as zero-momentum in continuum model 
theories, and separate quadratically with increasing wavevec- 
tor. In a IDES left and right going electrons cross the Fermi 
energy at a single point. In this figure the momentum of right- 
going (left-going) electrons is plotted relative to +kp (—ftp) 
where &f is the Fermi wavevector. 

Similarities and differences between graphene bilayers 
and IDES are most easily explained by temporarily ne- 
glecting the spin, and in the case of graphene also the 
additional valley degree of freedom. As illustrated in 
Fig. [1] in both cases the Fermi sea is point-like and there 
is a gap between occupied and empty free-particle states 
which grows with wavevector, linearly in the IDES case. 
These circumstances are known to support a mean-field 
broken symmetry state in which phase coherence is es- 
tablished between conduction and valence band states 
for arbitrarily weak repulsive interactions. In the case 
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of IDES, the broken symmetry state corresponds physi- 
cally to a charge density-wave (CDW) state, while in the 
case of bilayer graphene^ it corresponds to state in which 
charge is spontaneously transferred between layers. This 
mean-field theory prediction is famously incorrect in the 
IDES case, and the origin of the failure can be elegantly 
identified^ using a perturbative renormalization group 
(PRG) approach. We show below that when applied to 
bilayer graphcne, the same considerations lead to a dif- 
ferent conclusion. 

The reliability of the mean-field theory prediction^ of 
a weak-interaction instability in bilayer graphene can 
be systematically assessed using PRG4. We outline the 
main steps in the application of this analysis to bilayer 
graphene in the main text, pointing out essential differ- 
ences compared to the IDES case. Details are provided 
in the supplementary material. We assume short-range 
interactions* between electrons in the same (S) and dif- 
ferent (D) layers. 

The PRG analysis centers on the four point scattering 
function defined in terms of Feynman diagrams in Fig. [2] 
Since the Pauli exclusion principle implies that (in the 
spinless valleyless case) no pair of electrons can share the 
same 2D position unless they arc in opposite layers, in- 
tralaycr interactions cannot influence the particles; there 
is therefore only one type of interaction generated by the 
RG flow, interactions between electrons in opposite lay- 
ers with the renormalized coupling parameter Yd- The 
direct and exchange first order processes in Fig. [2] have 
the values Vd and respectively where Vd is the bare 
coupling parameter. 

The PRG analysis determines how Vd is renormalized 
in a RG procedure in which fast (high energy) degrees 
of freedom are integrated out and the fermion fields of 
the slow (low energy) degrees of freedom are rescaled to 
leave the free-particle action invariant. The effective in- 
teraction Td is altered by coupling between low and high 
energy degrees of freedom. At one loop level this inter- 
action is described^ by the three higher order diagrams 
labeled ZS, ZS', and BCS in Fig. [2] The internal loops in 
these diagrams are summed over the high-energy labels. 
In the case of IDES the ZS loop vanishes and the ZS' 
and BCS diagrams cancel, implying that the interaction 
strengths do not flow to large values and that neither 
the CDW repulsive interaction nor the BCS attractive 
interaction instabilities predicted by mean-field theory 
survive the quantum fluctuations they neglect. The key 
message of this paper is summarized by two observations 
about the properties of these one-loop diagrams in the 
bilayer graphene case; i) the particle-particle (BCS) and 
particle-hole (ZS, ZS') loops have the same logarithmic 



* We replace the bare Coulomb interactions by short-range 
momentum-independent interactional by evaluating them at typ- 
ical momentum transfers at the model's high-energy limit. We 
believe that this approximation is not serious because of screen- 
ing. 



divergences as in the IDES case in spite of the larger 
space dimension and ii) the ZS loop, which vanishes in 
the IDES case, is finite in the bilayer graphene case and 
the BCS loop vanishes instead. Both of these changes 
are due to a layer pseudospin triplet contribution to the 
single-particle Green's function as we explain below. The 
net result is that interactions flow to strong coupling even 
more strongly than in the mean-field approximation. The 
following paragraphs outline key steps in the calculations 
which support these conclusions. 




FIG. 2: Four-point scattering function for bilayer 
graphene. The 1st diagram is the renormalized interaction 
while the 2nd and 3rd ones are the direct and exchange bare 
interactions, respectively. The following diagrams are the one- 
loop diagrams labelled ZS, ZS' and BCS. The external and 
internal Green's function labels refer to layer in the case of 
graphene and to chirality in the case of IDES. 

An elementary calculation shows that the single- 
particle Matsubara Green's function corresponding to the 
Hamiltonian in Eq. {1} is 
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and hujk = £fc = Hk 2 /2m*. The pseudospin-singlet com- 
ponent of the Green's function which is diagonal 
in layer index, changes sign under frequency inversion 
whereas the triplet component ^t, which is off-diagonal, 
is invariant. 

The loop diagrams are evaluated by summing the prod- 
uct of two Green's functions (corresponding to the two 
arms of the Feynman diagram loops) over momentum 
and frequency The frequency sums are standard and 
yield (J3 = (fceT)- 1 ) 



^£^(-n) = T 



tanh(/% fc /2) _ .J_ 



(3K 



H 2 ^ T->0 
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where k is the momentum label shared by the Green's 
functions. Note that the singlet-triplet product sum van- 
ishes in the low-temperature limit in which we are inter- 
ested. Each loop diagram is multiplied by appropriate in- 
teraction constants (discussed below) and then integrated 
over high energy momentum labels up to the massive chi- 
ral fcrmion model's ultraviolet cutoff A: 



A/s<fc<A 



d 2 k 



tanh(/% fc /2) 
4^ 



(5) 



where vq = m* /2ttH 2 is the graphene bilayer density-of- 
states. Because lo^ oc k 2 , this integral grows logarithmi- 
cally when the high-energy cut-off is scaled down by a 
factor of s in the RG transformation, exactly like IDES. 
This rather surprising property of bilayer graphene is di- 
rectly related to its unusual band structure with Fermi 
points rather than Fermi lines and quadratic rather than 
linear dispersion. 

The key differences between bilayer graphene and the 
IDES appear upon identifying the coupling factors which 
are attached to the loop diagrams. The external legs in 
the scattering function Feynman diagrams (Fig. are 
labeled by layer index (T = top layer and B = bottom 
layer) in bilayer graphene. The corresponding labels for 
the IDES are chirality [R = right-going and L = left 
going); we call these interaction labels when we refer to 
the two cases generically. Since only opposite layer in- 
teractions are relevant, all scattering functions have two 
incoming particles with opposite layer labels and two out- 
going particles with opposite layer labels. In the ZS loop, 
at the upper vertex the incoming and outgoing T parti- 
cles induce a B particle-hole pair in the loop while the 
incoming and outgoing B particles at the lower vertex in- 
duce a T particle-hole pair. Since the particle-hole pairs 
must annihilate each other, there is a contribution only 
if the single-particle Green's function is off-diagonal in 
interaction labels. This loop can be thought of as screen- 
ing Vd; the sign of the screening contribution is oppo- 
site to normal, enhancing the bare interlayer interaction, 
because the polarization loop involves layer pscudospin 
triplet propagation. (See Eq. Q.) This contribution 
is absent in the IDES case because propagation is al- 
ways diagonal in interaction labels. The ZS' channel cor- 
responds to repeated interaction between a T particle 
and a B hole. This loop diagram involves only particle- 
propagation that is diagonal in interaction labels and its 
evaluation in the graphene bilayer case therefore closely 
follows the IDES calculation. This is the channel respon- 
sible for the IDES mean-field CDW instability in which 
coherence is established between R and L particles. In 
both graphene bilayer and IDES cases it has the effect 
of enhancing repulsive interactions. The BCS loop cor- 
responds to repeated interaction between the two incom- 
ing particles. In the IDES case the contribution from 
this loop which enhances attractive interactions, cancels 
the ZS' contribution, leading to marginal interactions 
and Luttinger liquid behavior. In the graphene bilayer 
case however, there is an additional contribution to the 



TABLE I: Summary of contrasting the contributions (in units 
of the related density-of-states) of the three one- loop diagrams 
in IDES and graphene bilayer cases 



diagrams 


ZS 


ZS' 


BCS 


one-loop 


IDES 





l? ln(s) 


-vf- ln(s) 





graphene bilayer 


±r£ ln( S ) 


i rj> ln( S ) 





Id ln(a) 



BCS loop contribution in which the incoming T and B 
particles both change interaction labels before the sec- 
ond interaction. This contribution is possible because of 
the triplet component of the particle propagation and, in 
light of Eq. (j4|) , gives a BCS loop contribution with a sign 
opposite to the conventional one. Summing both terms, 
it follows that the BCS loop contribution is absent in the 
graphene bilayer case. These results arc summarized in 
Tabic U and imply that at one loop level 



1 - Vdvq ht(s) 



(6) 



The interaction strength diverges when Vd vq = 1/ ln(s), 
at halft the mean-field theory critical interaction 
strength. Taking guidance from the mean- field theory 5 , 
the strong coupling state is likely a pscudospin fcrromag- 
net which has an energy gap and spontaneous charge 
transfer between layers. 

A number of real-world complications which have to 
be recognized in assessing the experimental implications 
of these results. First of all, electrons in real graphene 
bilayers carry spin and valley as well as layer pscudospin 
degrees of freedom. This substantially complicates the 
PRG analysis since many different types of interactions 
are generated by the RG flow. When one spin degree 
of freedom is considered three types of interactions have 
to be recognized, Tg, Tp, and Fx- couples electrons 
with the same flavor and Td electrons with different fla- 
vors, while Tx is the exchange counterpart of Tjj. The 
bare value of the exchange part of (the direct part 
of Tx) is of course zero since the Coulomb interaction 
is flavor independent, but higher order contributions are 
nonzero if triplet electron propagation is allowed. Corre- 
spondingly, when both spin and valley degrees of free- 
dom are acknowledged the interaction parameters are 

TsSD, TsDS, TsDD, TdSS, ^DSD, ^DDS, ^DDD, ^XSD, 

Txds, and Txdd (see supplementary material) . The 
labels refer from left to right to sublattice pseudospin, 
real spin, and valley degrees of freedom. Tsss is ab- 
sent due to Pauli exclusion principle while Txss is just 



' Mean-field theory is equivalent to a single-loop PRG calculation 
in which only one particle-hole channel is retained. For the draw- 
ing conventions of Fig. [2] the susceptibility which diverges at the 
pscudospin ferromagnet phase boundary is obtained by closing 
the scattering function with r z vertices at top and bottom, so 
the appropriate particle-hole channel is the ZS channel not the 
ZS' channel. 
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—Tdss due to the fermionic antisymmetry between out- 
going particles. The one loop flow of these ten interaction 
parameters is described in the supplementary material. 
We find that the renormalized interactions diverge near 
Vd — 0.6/ln(s). The instability tendency is therefore 
somewhat enhanced by the spin and the valley degrees 
of freedom. 

Next we must recognize that the massive chiral fcrmion 
model applies only between trigonal-warping and inter- 
layer hopping energy scales. Appropriate values for ln(s) 
in the RG flows are therefore at most around 2.3 (see sup- 
plementary material) in bilayer graphene, compared to 
the unlimited values of the chiral-fermion model. When 
combined with estimates of the bare interaction strengths 
(see supplementary material), this limit on interaction 
flow suggests that spontaneous gaps are likely in sus- 
pended bilayer graphene samples, and much less likely 
for graphene bilayers on substrates. 

Gaps do appear in bilayer graphene even when 
electron-electron interactions are neglected, provided 
that an external potential difference V is applied be- 
tween the layers. The potential adds a single-particle 
term — Vt z /2 to the single-particle Hamiltonian, breaks 



inversion symmetry, and transfers charg e 10 ' 11 ^ 2 ' 13 be- 
tween layers. This interesting property is in fact the 
basis of one strategy currently being explored in the ef- 
fort to make useful electronic device a 14 i 15 out of graphene 
2DES's. Even if gaps do not appear spontaneously in real 
bilayer graphene samples, it is clear from the present 
work that interesting many-body physics beyond that 
captured by commonly used electronic-structure-theory 
approximations (LDA or GGA approximations for exam- 
ple), must play at least a quantitative role in determining 
gap grown with V. As graphene bilayer sample quality 
improves, we expect that it will be possible to explore 
this physics experimentally with ARPES, tunneling, and 
transport probes. 
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In this supplementary information, we summarize the Green's function technique, detail the diagrammatic pertur- 
bation theory calculations, and explain the many interaction channels generated by the RG flows when spin and valley 
degrees of freedom are included. We also discuss the bare interactions and the important trigonal-warping effect. 

I. GREEN'S FUNCTIONS AND FREQUENCY SUMS 

The low-energy effective Hamiltonian of a graphene bilayer is a momentum q-dependent 2x2 matrix of the following 
form: 

34? = —B ■ t , (1) 

where r is the Pauli matrix vector, |B| = £ q = hcj q = h 2 q 2 /2m*, and the orientation angle of B is twice the orientation 
angle of q (</> q ). (This property is often expressed by saying that the model has layer pseudospin chiralityi^ J=2.) 
The pseudospin degree of freedom upon which the operator t acts represents the layer in which an electron resides. 
Eigenstates of Jif are coherent linear combinations of amplitudes in both layers. The finite temperature Green's 
function at Matsubara fermion frequency iuj n and wavevector q is given by 

Sf (iw n ) = [iw„ - Jif/h]- 1 = %{iu n ) + Sf t («w„) ti t (2) 

where 

%^n) = \ {— ± . \ ) , (3) 

n = —B/\B\, the "s" and "t" subscripts denote singlet and triplet contributions respectively, and we have chosen 
jit = to address the properties of an electrically neutral bilayer. Note that Sf s (— ioj n ) = — & e (iw n ) whereas §^(— iu) n ) = 
^t(w n ). Using the property that n = — (cos(20 q ), sin(20 q ), 0) we obtain the following explicit form for the 2x2 
Green's function matrix: 

The off diagonal triplet component captures processes in which free electrons propagate from one layer to the other. 
Evaluation of the loop diagrams which appear in the RG calculation described in the main text requires frequency 
sums to be performed for products of two Green's functions. These are given by 

-Lycg? t (iu n ) = T tanh (^/ 2 ) _^ =p— i— , *Y9.(iu n )& t (iu> n )^>0. (5) 

P h ^ T "° ^ q P % zt 

II. SPINS PSEUDOSPINS AND DISTINCT INTERACTION PARAMETERS 



In the low-energy continuum model of bilayer graphene electrons carry spin, and both layer and valley pseudospin 
labels. In a scattering event, both the two incoming and two outgoing particles can therefore have one of eight 
labels and the general scattering function therefore has 8 4 possible low-energy long- wavelength values. The number 
of distinct coupling constants in the RG flow equations is much smaller, however, because many values are zero and 
others are related to each other by symmetry. One simplification is that interactions conserve spin, and both layer 
and valley pseudospin, at each vertex. Interactions are however dependent on whether the interacting particles are 
in the same (S) or in different (D) layers. The internal loops in the perturbative RG calculation contain two fermion 
propagator (Green's function) lines. These propagators conserve both spin and valley pseudospin, but as we have 
seen above, not the layer pseudospin. It is clear then that the incoming and outgoing total spin must be preserved for 
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FIG. 1: Electron-electron scattering processes for a system with one pseudospin-1/2 degree of freedom. 

real spin and for the valley pscudospin, but the layer pseudospin case requires a more elaborate consideration. From 
Eq. ([4| we sec that a phase factor e ^ is gained when the propagator transfers electrons between layer index with 
the + for top to bottom evolution and the — for bottom to top. Unless these transfers enter an equal number of times 
in each direction, the integrand in a Feynman diagram will contain a net phase factor related to chirality and vanish 
under momentum integration. The total layer pseudospin is therefore also conserved in collisions. 

In identifying distinct coupling constants, we start with the simplest case in which the valley and spin labels are 
absent. There arc then three possibilities, as illustrated in Fig. [TJ When the two incoming pscudospins are parallel 
(r,5 in Fig. [TJ) , the outgoing pscudospins must also be parallel. Because of Fermi statistics interchanging the outgoing 
lines in T$ changes the diagrams's sign. Since the diagram is invariant under this operation, it must vanish. The 
second possibility is opposite incoming pscudospins, which requires opposite outgoing pscudospins in one of the two 
configurations labelled by Tp and Fx in Fig.[TJ In this case Fermi statistics implies that To = —Tx- It follows that 
the only distinct interaction parameter is Tp. 

If more than one pseudospin is present, we have to recognize more separate interacting processes. For example, 
for systems with two relevant pseudospins, the interaction parameters can be labeled in the same way as in Fig. [T] 
but by doublets which account for the different pseudospins separately Two pscudospin interactions might include 
^sd, r_DS- Tqd and Txd (see Fig. [2]) for example. In models for which propagators and interactions preserve all 
pseudospin labels, we would have Txd — since all pscudospin flavors are preserved along each fermion line. For 
graphene bilaycrs, however, we must keep Txd ^ because the layer pseudospin has triplet propagation. Following 
this line of argument, we conclude that in graphene bilayers, with its three different pseudospins, there are ten distinct 
non-zero interaction parameters: T S sd, T S ds, T S dd, T DS s, T DS d, T DDSl T DDDl T X sd, T X ds and T X dd, where 
the first label refers to layer pseudospin, and the following labels to real spin and valley. 




FIG. 2: Distinct interaction channels for systems with two pseudospin-1/2 degrees of freedom. In this figure the 
first spin is denoted by | f) or j) while the second by | — >) or | <— ). The one-dimensional electron gas system can be viewed 
as being in this class if the chirality index is regarded as a pseudospin. 



III. PRG ANALYSIS FOR SPINLESS AND VALLEYLESS GRAPHENE BILAYERS 

The one-loop correction to the interlayer interaction Tp can be calculated using the following ZS, ZS' and BCS 
diagrams^. Since Tp (with bare value Vd) is the only interaction parameter as we explained previously, the propagator 
labels of the internal loops must be distributed accordingly. We label the external legs in the scattering diagrams by 
their layer indices (T = top layer and B = bottom layer). The corresponding labels for the IDES are chirality (R = 
right-going and L = left going). 
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FIG. 3: ZS loop correction in the one-loop perturbative RG calculation. 



In the ZS loop shown in Fig. [31 at the upper vertex the incoming and the outgoing T particles induce a B particle- 
hole pair in the loop while the incoming and outgoing B particles at the lower vertex induce a T particle-hole pair. 
The corresponding labels in the IDES case arc L, R for left and right chirality. This correction survives in the bilayer 
case only because the single-particle Green's function has a triplet contribution [see Eq. ((5])] which is off-diagonal in 
layer index. The ZS contribution is absent in the IDES case^ because propagation is always diagonal in interaction 
labels. Here we find 



^2 

^ZS 
D 



~ l.^(Wn)-r D J — --T D v ln(s), (6) 

where vq = m* /2itTi 2 is the graphene bilayer density-of-states (per spin and valley) and the integral is carried out in 
the momentum shell A/s < q < A. In the zero temperature limit 



A/s<g<A 



l q tanh(/?£j/2) 1 

' 7^0 ln ( S ) • ( 7 ) 



(2tt) 2 4£ q t^o 2 



The ZS' channel shown in Fig. 3] corresponds to repeated interaction between a T particle and a B hole. This 
loop diagram involves only particle-propagation that is diagonal in interaction labels; its evaluation for the IDES 
and graphene bilayers correspond closely. This is the channel responsible for the IDES mean-field CDW instability^ 
in which coherence is established between R and L particles^. In both cases it has the effect of enhancing repulsive 
interactions. We find 




FIG. 4: ZS' loop correction in the one-loop perturbative RG. 



r zs' r l f d2 Q v-«2, • n P 2 f d2 Q tanhQgg/2) 1 2 

Td = ~0?J pL^^^) = r fl / — ^ — = -r D v Hs). (8) 

The BCS loop corresponds to repeated interaction between the two incoming particles. In the IDES case the 
contribution from this loop (see Fig. [5]) cancels the ZS' contributions, leading to marginal interactions and Luttingcr 
liquid behavior. This same kind of BCS correction for graphene bilayer reads 
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FIG. 5: BCS (particle-particle) loop correction for singlet propagation in the one-loop perturbative RG calcu- 
lation. 
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/t2 



2ph 2 J (2tt) 

i(-r D ) 2 f d 2 q 



2 ^ 2 7 (2^) 2 ^ 
-T 2 D v ln(s) . 



,(q,iuj n )%(-q, -iu n ) 



(9) 



In the graphene bilayer case, however, there is an additional contribution (see Fig. [5]) to the BCS loop contribution 
in which the incoming T and B particles both change layer labels before the second interaction. This contribution 
is possible because of the triplet layer pseudospin propagation and, in light of Eq. (O, gives a BCS contribution of 
opposite sign to the normal contribution: 
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FIG. 6: BCS(particle-particle) loop correction for triplet propagation in the one-loop perturbative RG calcula- 
tion. 
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It follows that the BCS loop contribution is absent in the graphene bilayer case because 



r 



BCS 



r 



BCSi 



r 



BCS 2 



0. 



D " 1 D "T" x D 

Therefore, at one-loop level, the renormalization of intcrlayer interaction is 



^one — loop 
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rg* + if»' +r^« = iVoin0 5 ): 



(10) 

(ii) 

(12) 
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which leads to the RG flow equation 



dT 



D 



vo dln(s) 

Combined with the bare interlayer interaction Vd and integrating the flow equation we obtain that 

Vd 



D 



(13) 



(14) 



1 - VD^ohi(s) 

which diverges at the point Vo fo — 1/ hi(s). 

For the Fcynman diagrams drawing conventions we have chosen, the interaction correction to the layer pscudospin 
response function \ zz , which diverges at the pseudospin ferromagnet phase boundary is obtained by closing the 
scattering function with a t z vertex at top and bottom. The t z operator measures the charge difference between T 
and B layers. Because it is an effective single-particle theory, fermion mean-field theory^ corresponds to response 
function diagrams with at most a single particle-hole pair. It follows that mean-field theory is equivalent to a single- 
loop PRG calculation in which the BCS and ZS' channels are neglected and only the ZS channels is retained. In 
mean-field theory^ the ideal graphcne bilayer has an instability to a state in which charge is spontaneously transferred 
between the layers which is signalled by the divergence of Xzz- The PRG analysis demonstrates that the mean- field 
theory instability is enhanced by reinforcing ZS' channel contribution. 



IV. PRG ANALYSIS FOR SPINFUL AND VALLEYFUL GRAPHENE BILAYERS 



As discussed at the end of section II, there are ten interaction parameters for the spinful and valleyful case. The 
one-loop flow equations are derived in the same way as in the spinless valleyless case, except for the necessity of 
keeping track of the many-possible configurations of the end labels on the loop propagators. We find that 
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dT SDS 
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vq rfln(s) 

dToDD 

i/q dln(s) 
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vq rfln(s) 

dTxps 
vq dln(s) 

dTxDD 

vq dln(s) 
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- Fdss(Fddd — Tsdd) — (Tdsd — ^ssd)(^dds — Tsds) 



(Txsd — FssdY 



(T 



XDS 



SDS) 



r 

ir 



1 



1 



DSS + o^DSD — FsSd) + -^(TdDS - F,SDs) + DDD ~ ^SDd) 



1, 



XSD 



XDS 



DSD 

— r 2 

2 1 DDS 



DSS^DSD 



-r 



SSD) 



(XDDD — ^SDd)(XdDS — ^sds) — 2 O^DSD 



FdSs(TdDS — ^sds) + (X DDD — ^SDd)(^DSD — ^SSd) — ^i^DDS 
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DDS 
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SSD) 



DDD 
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FxDs) 2 
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FdSS^XSD — FxDsFxDD — r (TxSD — FsSd) 2 — ^(^ X SD + ^DSd) 2 



Fx 



DS 
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1 2 1 

TdSS^XDD — ^XSD^XDS — 2^ XDD ~ ^SDD) ~ -^{^XDD 



DDS) 



DDD) 



XDD 



(15) 



The only fixed point that we have identified is the non-interacting one. These ten coupled flow equations can be 
integrated numerically starting from bare interactions. In order to represent the property that same layer interactions 
will be slightly stronger than different layer interactions we set the bare interactions values to 1.1,0.9 and for 

Vssd (V S ds,Vsdd) ,Vq V D ss {Vdsd,V DDSi V ddd ) and v V X sd (V X ds,V X dd), respectively. (The motivation for 
this choice is explained in the next section.) We find that the interaction parameters flow away from the non-interacting 
fixed point and diverge at a finite value of s as illustrated in Fig. [7] The instability criterion implied by this one-loop 
PRG calculation is Vo v$ ~ 0.6/ln(s). The instability tendency is therefore enhanced by the spin and valley degrees 
of freedom since the criterion was Vo vq — l/ln(s) for the spinless and valleyless case. 
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FIG. 7: Renormalized interaction flow. This illustration plots the inverse interaction strength (fore e ff) _1 versus the scaling 
parameter ln(s)/e e ff. e e fl is the effective dieletric constant of the graphene bilayer and V = IVacuum/seff. Interlayer interaction 
parameters Tdss (green) and Fxsd, Fxds, Vxdd (cyan) flow to large values most quickly. According to this estimate the 
normal state becomes unstable for Vd larger than ~ 0.6/ m(s). 



V. BARE INTERACTIONS AND THE TRIGONAL- WARPING EFFECT 



The conclusions which can be drawn from the PRG calculation presented here are sensitive to the upper and lower 
momentum and energy cutoffs, which limit the applicability of the massive chiral fermion model for bilayer graphene, 
and to the strength of bare electron-electron scattering amplitudes. Below we estimate the range of s over which the 
RG flows discussed above apply, and the strength of the bare interaction Vd- We caution that, given the nature of 
the PRG calculations, the estimates presented below should be regarded as qualitative. 

In practice the upper cutoff is the interlayer hopping energy 71 ~ 0.4 eV; at higher energies it is essential to account 
for two sublattice sites in each layer. We have in addition ignored the trigonal-warping part in the full Hamiltonian, 
due to a direct hopping process between the low-energy sites which has energy scaled 73 ~ 0.3 cV. Inserting the 
expression^ for the effective mass of the massive chiral Fermion model we find that the model we have studied is 
appropriate for 

73 H 2 q 2 H 2 q 2 v 2 

hv F q— < - — - = -<7i (16) 

7o 2 m* 71 

where i>f ~ 10 8 cm/s is the Fermi velocity near the Dirac point in the single-layer-graphene continuum model, 
and 70 ~ 3 eV is the intralayer near neighbor hopping energy. It follows that the high energy momentum cutoff 
5h = 7i/ft«F an d that the low energy momentum cutoff = (7371/70)/^^, which gives the maximum value of the 
scaling parameter ln(s). Using accepted values for the hopping parameters^, it follows that the scaling relations we 
derive should apply approximately over a wavevector range corresponding to ln(s) max = \n(qn / q-^) ~ 111(70/73) ~ 2.3. 

We estimate the strength of the bare scattering amplitudes by evaluating the 2D Coulomb scattering potential at 
the cut-off wavevector feu : 

to* 27re 2 a ec 

where a cc — e 2 /hvp ~ 2.2 is graphene's fine structure constant. The value used for Vg in the RG flows is motivated 
by this estimate. The value used for Vd is reduced by a factor of exp(— k^d) compared to Vs to account for the layer 
separation d = 3.35 A. 

According to these estimates the bare value of VqVd exceeds the stability limit of ~ 0.6/ ln(s) max ~ 0.25 by 
approximately a factor of four. The above estimates are for the case of a graphene bilayer in vacuum. For graphene 
layers on the surface of a substrate with dielectric constant e, interactions are expected to be reduced by a factor of 
~ (e + l)/2. In the case of SiO-2 substrates e ~ 4 and the interaction strength exceeds the stability limit by a much 
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narrower margin. Wc expect that additional screening effects from graphene a orbitals, which arc normally neglected 
in continuum model calculations, will reduce interaction strengths at wavevectors near fcjj somewhat and favor stable 
bilayers. 
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